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Abstract: We present the O(aj) virtual QCD corrections to unlike-quark qq — > q'q' 
and like-quark scattering qq — > qq due to the interference of one-loop amplitudes with 
one-loop amplitudes. The structure of the infrared divergences agrees with that pre- 
dicted by Catani. The results are expressed in an analytic form so that the relevant 
expressions for crossed scattering processes can be obtained in a straightforward 
manner. The one-loop contributions presented here, together with the interference 
of tree with two-loop amplitudes given recently, complete the 2^2 virtual QCD 
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1. Introduction 



The increasing precision of high energy scattering experiments demonstrates the 
need for very precise theoretical calculations. The accuracy of existing next-to- 
leading-order (NLO) predictions may be improved by including the next highest term 
in the perturbation series. Such next-to-next-to-leading-order (NNLO) estimates 
will improve the theoretical precision in two ways. First, the renormalisation scale 
dependence will decrease (from about 10% for NLO predictions for central production 
of a jet with transverse energy E T ~ 100 GeV at the Tevatron to approximately 
1 — 2% at NNLO). Second, better matching of the parton level theoretical and hadron 
level experimental jet algorithms will be possible since at NNLO, three partons can 
combine to form the jet rather than two. 

Recent progress towards two- loop integrals [l]-[8] has turned the calculation of 
NNLO virtual corrections for massless 2—^2 processes into a viable task. Bern, 
Dixon and Kosower [9] were the first to address such scattering processes and pro- 
vided analytic expressions for the maximal-helicity-violating two-loop amplitude for 
99 — * 99- More recently, Bern, Dixon and Ghinculov [10] completed the two- 
loop calculation of physical 2 2 scattering amplitudes for the QED processes 
e + e~ — > and e + e~ — > e~e + . 

Subsequently, we have derived analytical expressions for the O(aj) two-loop 
contribution to unlike quark scattering qq — > q'q' [11], and like quark scattering qq — > 
qq [12], as well as the crossed and time reversed processes, in the limit where the quark 
mass can be neglected. To complete the calculation of the NNLO virtual corrections, 
the interference terms of one-loop amplitudes with one-loop amplitudes need to be 
included. It is the purpose of this paper to present analytical expressions for these 
contributions using conventional dimensional regular isat ion (space-time dimension 
D = 4 — 2e) , renormalised with the MS scheme. 

Our paper is organised as follows. We establish our notation in Section 2 and 
briefly describe our method in Section 3. In Sections 4.1 and 4.2, we provide analytic 
expressions of the one-loop contributions at NNLO for the unlike and like-quark 
scattering respectively, obtained by direct evaluation of the Feynman diagrams. Our 
results are expressed in terms of two master integrals, the one-loop box graph in 6— 2e 
dimensions and the one-loop bubble graph in 4 — 2e dimensions. A clear separation 
of the infrared poles is apparent by inspection of the results. We find complete 
agreement between the infrared structure obtained by our explicit calculation and 
that anticipated by Catani [13]. Analytic expansions in e for all kinematic regions 
are straightforward to derive by inserting the expansions of the master integrals in 
the appropriate region. In Section 5 we summarize our results. 
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2. Notation 



We follow the notation of Refs. [11, 12] as closely as possible and consider the unlike- 
quark scattering process 



q(pi) + q(p2) + q'(ps) + q'{pi) -»• o, 

and the like-quark scattering process 

q(pi) + q{p2) + q(ps) + q{pi) -> o, 



(2.1) 



(2.2) 



where particles are incoming and carry light-like momenta (shown in parentheses). 
Their total momentum is conserved, satisfying 



P?+P%+P%+Pl = 0, 
and the associated Mandelstam variables are given by 

s = (pi + P2) 2 , t = (p 2 + p 3 ) 2 , u =(p 1 + p 3 ) 2 . 



(2.3) 



We use conventional dimensional regularisation and treat the external quark states 
in D space-time dimensions and renormalise the ultraviolet divergences in the MS 
scheme. The bare coupling a is related to the running coupling a s = a s (fi 2 ), at 
renormalisation scale ji, by 



a S e = a s 



A) ( ot 



2tt 



+ 



'P 2 o Pi 



2e V2tt/ 



+ 0{al) 



where 



(47r) e e" 



■e-y 



7 = 0.5772 



Euler constant, 



(2.4) 



(2.5) 



is the typical phase-space volume factor in D = 4 — 2e dimensions. As usual, the 
first two coefficients of the QCD beta function, (3 and f3\ for N F (massless) quark 
flavours are 



A) 



IIC A -AT R N F 



Pi 



17C\ - IOCaTr N f - 6C F T R N F 
6 



where iV is the number of colours, and 

N 2 -l 



C A = N, Tj 



R 



2N 1 2 

The renormalised amplitude for the unlike-quark process is given by 



\M) 



unlike 



Ana. 



|^ (0) > + (^) l^ (1) >+(^) 2 |^ (2) ) + 0(^) 



(2.6) 



(2.7) 



(2.8) 
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with |.M^) representing the i-loop amplitude in colour-space. For the like-quarks 
we have the related expression 



\M)uke = 4vra 5 



|A<< >) - \M i0) ) ) + g ) ( \M^) - \M (1) ) 

+ (S 2 {\^)-\M^)) + { 



«2) 



(2.9) 



u\ 

Here \A4 ) describes the t-channel graphs which can be obtained from the s-channel 
diagrams by exchanging the roles of particles 2 and 4 

\M {€) ) = \M (i) )(2 <->4). (2.10) 

Both \M^) and \M^) are renormalisation scale and renormalisation scheme depen- 
dent. 

In squaring the amplitudes and summing over colours and spins we find two 
types of terms, 

• the self-interference of the graphs in a single channel, described by the function 
A(s, t, u) for the s-channel and A(t, s, u) for the t-channel, and 

• the interference of the s-channel graphs with the t-channel graphs, described 
by the function B(s,t,u). 

Thus, for distinct quark scattering we have 

(M\M) unHke = ]T \M(q + q - q' + q')\ 2 = A(s, t, u), (2.11) 

while for identical quarks 

(M\M) Hke = J2\M(q + q^q + q)\ 2 

= A(s, t, u) + A(t, s, u) + B(s, t, u). (2.12) 

Similarly, for the crossed and time-reversed processes we obtain 

]T \M(q + q' - q + q')\ 2 = A(u, t, s) (2.13) 

J2\M(q + q'^q + q')\ 2 = A(t, s, u) (2.14) 

£ \ M (q + 4 - q + q) ? = A(u, t, s) (2.15) 

\M(q+ q^q + q)\ 2 = A(u,t, s) + A(t,u, s) + B(u,t, s). (2.16) 

The function A can be expanded perturbatively to yield 



A(s, t, u) = 16ir 2 a 2 



A%s, t, u) + A%s, t, u) + (g) 2 A s (s, t, u) + O(al) 



(2.17) 
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where 



A\s,t,u) = (MW\M<®) = 2(N 2 - 1) ( - e) , (2.18) 



A\s,t,u) = ((M {0) \M {1) ) + (M (1) \M (0) )) , (2.19) 
A 8 (s,t,u) = ({M {1) \M {1) ) + (M {0) \M (2) ) + (M^\M {0) )) . (2.20) 

In the same manner 



B(s,t,u) = 16n 2 a 2 



B 4 (s, t, u) + (g) B"(s, t, u) + (g) 2 £ 8 ( S , f , u) + O(al) 



(2.21) 

where, in terms of the amplitudes, we have 
B%s,t,u) = - (^(M {0) \M^) + (M^\M {0) )^J 

B%s,t,u) = - ((M (1) \M (0) ) + (M (0) \M {1) ) + (M (0) \M {1) ) + (M W \M {0) ) 

(2.23) 

B 8 (s,t,u) = - ((M (1) \M (1) ) + (M (1) \M (1) ) 

+ (M (0) \M^) + (7W^|7W (0) > + (M^\M {2) ) + (7W (2) |^(°))) . 

(2.24) 

Expressions for A 6 and £> 6 , valid in dimensional regularisation, are given in Ref. [14]. 
Analytical formulae for the two-loop contribution to A 8 

(M i0) \M {2) ) + (M^\M^) 

are given in Ref. [11], with analogous expressions for the two-loop contribution to B 8 

(M i0) \M^) + (M^\M m ) + (M^\M {2) ) + (M (2) \M^) 

given in Ref. [12]. 

We now concentrate on the contributions to both A 8 and B 8 due to the interfer- 
ence of one-loop amplitudes with one-loop amplitudes, namely 

A 8(lxl \s,t,u) = (M {1) \M {1) ), (2.25) 

and 

B 8{lxl \s,t,u) = - ((M (1) \M (1) ) + (M W \M {1) )^ . (2.26) 
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Even though they are somewhat simpler to evaluate than the two loop graphs, they 
form a vital part of the NNLO virtual corrections and we present them here for 
completeness. One-loop helicity amplitudes for the 2 — > 2 quark scattering processes 
were given in Ref. [15] as truncated expansions in e including their finite part. How- 
ever, this is only sufficient to obtain the pole structure of v 4. 8 ( lxl ) and £> 8<lxl ) up 
to 1/e 2 . To determine the 1/e and finite parts requires knowledge of the one-loop 
amplitude through to 0(e 2 ). 



3. Method 

We organise our calculation as follows. First the one-loop Feynman diagrams are 
generated using QGRAF [16]. The emerging tensor integrals are associated with scalar 
integrals in higher dimension and with higher powers of propagators [6, 17]. Sys- 
tematic application of the integration-by-parts identities [18] is sufficient to reduce 
these higher-dimension, higher-power integrals to master integrals in D = 4 — 2e. We 
can express all the integrals of the one- loop amplitudes in terms of just two master 
integrals, the scalar bubble graph 

Bub( S )= -Q-(s) , 

and the one-loop scalar box graph 



Box(s,t) = (s, t) 



The above choice of master integrals is not unique. We prefer to replace the one-loop 
box in D = 4 — 2e by the finite one- loop box in D = 6 — 2e, Box 6 . This leads to a 
natural separation of the infrared poles and the finite part of the one-loop amplitudes. 
As a last step we multiply together the one-loop amplitudes and perform the colour 
and Dirac traces. 

4. Results 

In this section we give explicit formulae for both A 8( - lxl \s, t, u) and B 8 ( - lxl ^(s, t, u), 
in terms of the finite Box 6 and the 1/e divergent Bub master integrals. The finite 
parts depend only on Box 6 and differences of the Bub master integrals. 

4.1 Unlike quarks 

In the unlike-quark case we obtain, 

A 8{lxl \s,t,u) = [\in t + T r + T g \ 2 + (N 2 - l) \in nt \ 2 ] (Mo\M ) 



+ 2 Re 



(IKt + F r + + (N 2 - l) XTl{ t T 2 
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+ (n 2 - l) 



iV 4 — 3iV 2 + 3 _ , , iV 2 + 3 

-^ 3 (s, t, u) + F z {s, u, t) 



+ 



N 2 
N 2 -3 
N 2 



N 2 



[JF 4 (s, t, u) + F±{s, M, t)] 



(4.1) 



where the infrared poles present in the one-loop amplitude proportional to the tree- 
level matrix elements are given by 



in 



e(2 + e) 
2 



nt 



e(2 + e) 

which diverge as 1/e 2 and 1/e respectively. Both 

1 3(1 - e) 



1 2 (N 2 — 2) 

-Bub(s) - -Bub(u) - L__J B ub(0 

l B ub(u) - ^Bub(t) 



(4.2) 
(4.3) 



and 



e 3 - 2e 
e[iV 2 (ll + 2e) + 9-4e 2 



Bub(s) , 



Bub(s), 



2(2 + e)(3-2e)iV 

are finite terms multiplying the tree-level matrix elements. The functions 

iV 2 -l r 



?1 



2N 



(N 2 -2)f(s,t,u) + 2f(s,u,t) 



and 



N 2 — 1 
2iV 



f(s,t,u) - f(s,u,t) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



are finite and multiplied by the infrared poles of the conjugated one-loop amplitude, 
with 



f(s,t,u) 



3s 2 + 7u 2 + 9t 2 A u 2 + t 2 + 2s 2 5u + It 
4^; — h e 



(2 + e)s 2 

x 6t 2 + 2m 2 - 3es 2 „ ,, . 
+u(l-2e) = Box 6 (s,t). 



[Bub(t) - Bub(s) 



(4.8) 



Finally the square of the finite part of the one-loop amplitude is fixed by the finite 
functions JF 3 and JF 4 , 



F 3 (s,t,u) = |Box 6 (s,t) 
+ 2 Re 



2^~ 



+ |Bub(t) - Bub(s 



Bub(t) - Bub(s)J Box 6 (s,t) 
5t 2 - 2tu + 2m 2 



2m 3 - tu 2 + 8t 2 u - t 3 



2s 2 



+ 0(e), 



(4.9) 
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and 



F 4 (s,t,u) = 2Re{Box 6t (s,t)Box 6 (s,M)} 



tu(t 2 + u 2 



+2Re{[Bub(w) - Bub(s)] t Box 6 (s,t)} 



u(7t 2 - 2tu + 3u 2 ) 



+2Re{[Bub(w) - Bub(s)] f [Bub(f) - Bub(s)]} 



2s 2 

3(t 2 -tu + u 2 ) 
2^2 



+ 0(e). 
(4.10) 



In the latter expressions, we have discarded contributions of 0(e). 

After explicit series expansion in e, the infrared singular terms TlZ t and 11Z nt 
reproduce the pole structure obtained by expanding 



^7 



111 



t,C 



r 1-e 



1 3 
~n + 



3 e 7 



in 



nt,C 



1 



+ 



1 

N 
1 

TV 



.fL. 
s 

u 



2 
N 

1 

N 



u 
t 



(N 2 - 2) 
N 



' t 



(4.11) 

r(l-e) Ve 2 ' 2tJ [N \ u) N \ t J \' 

which is the singular structure obtained by straightforward application of the for- 
malism of [13]. To rewrite Eq. (4.1) directly in terms of 1TZt,c and XTZ n t,c rather 
than XTZ t and XlZ nt requires the finite difference to be evaluated through to 0(e 2 ). 

Equation (4.1) is valid in all kinematic regions. Series expansions in e in a 
particular region can be easily obtained by inserting the appropriate expansions of 
the master integrals. In this equation, the finite functions are multiplied by poles in 
e, so they must be expanded through to 0(e 2 ). In the physical region u < 0, t < 0, 
Box 6 (w,t) has no imaginary part and is given by [10] 



Box 6 (M) = 



x 



e^T(l + e)r(l-e) 2 / 
2sY (2 - 2e) \ s 

( 1 7T 2 

+2e I Li 3 (x) - L x Li 2 (x) - -L 3 X - —L, 



2e 2 (Li 4 (:r) + Lj,Li 3 (x) 



-((L x -L y ) 2 + n 2 ) 



7T 4 \ 



7T 2 IT 2 

—± L l ~ Y LxL y ~~ is J + (m ~ 



(4.13) 



where x = —t/s, L x = log(x) and L y = log(l — x) and the polylogarithms Li n (z) are 
defined by, 



lA n (z) = / — Li n _i(t) for n = 2, 3, 4 
Jo t 

, . f z dt , . 

Li 2 z = -/ -log 1-0- 

JO t 



(4.14) 
(4.15) 
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Analytic continuation to other kinematic regions is obtained using the inversion 
formulae for the arguments of the polylogarithms (see for example [6]) when x > 1, 

r 2 



Li 2 (x + iO) = -Li 2 



1 7T' 

-log 2 (x) + y + inlog(x) 



7T 



?7T 



Li 3 (a; + iO) = Li 3 - - - \og 6 (x) + — log(x) + — log 2 (a;) 
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U 4 ( X + iO) = -Li 4 (I) - 1 log 4 (:r) + y log 2 (x) + ^ + y log 3 (x) (4.16) 
Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by, 



Bub(s) = 



r (2 - 2e) e \ s 



(4.17) 



and can be easily expanded in e in all kinematic regions. Note that the complex 
conjugate of the master integrals is also required in Eq. (4.1). 

4.2 Like quarks 

For the like-quark contribution we find a similar expression, 
-2 Re | (TR t + T r + y 9 ) ] {TK t + F r + T g ) (M \M ) 

+ (N 2 - l)(TJZ nt -Tr~ TgfxTlnt (M \M ) 

(m t +T r + T g ) V x + (n 2 - \)(m nt -T r - T g ) ] r' 2 

fl(s,t,u)f 4 (t,s,u) + (s <-> f) 



+ 



(s <-> t) 



N 2 -l 

AT - 



+ 



N A — N 2 — 1 
2iV^ 
iV 4 - 2N 2 - 1 



27V 2 
3N 2 + 1 



2iV 2 



fl(s,t,u)f 4 (t,s,u) 



The infrared singular functions are given by 

2 



in 



e(2 + e) 
2 



1 1 (N 2 + l) 

-Bub(s) + -Bub(f) - ^pBub( M ) 

iV 2 — 1 1 1 

___Bub( S ) - yBub(t) + -Buh(u) 



e(2 + e 

which diverge as 1/e 2 . The finite renormalisation term is 



(4.18) 

(4.19) 
(4.20) 



— A) 



1 3(1-6) 
e 3-2e 



Bub(t)) , 



(4.21) 
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while the remaining finite contribution multiplying tree-level is given by 

— e(N 2 (ll + 2e) + 9 -4e 2 ). 



Fa 



-Bub(t). 



(4.22) 



9 2(2 + e)(3-2e)iV 

Once again, the finite part of the crossed one loop amplitude multiplying the infrared 
divergent terms of the one loop amplitude generates finite functions 

N 2 -l 



2N 2 



(N 2 -2)f 1 (s,t,u) + 2f 2 (s,t, 



u 



and 



where 



To 



N 2 -l 
2N 2 



[fi(s,t,u) - f 2 (s,t,u)] , 



(4.23) 
(4.24) 



fi(s,t,u) 



2u 

it 1 



1 - 2e) [t 2 + u 2 - 2e(t 2 + s 2 ) + e 2 s 2 J Box 6 (s, t) 

2u{2u -t) + e(u 2 -tu- At 2 ) + e 2 (ts - 4m 2 ) 
[Bub(t) - Bub(s)] , 



st{2 + e) 



+e 3 us + e ts 



(4.25) 



and 



f 2 (s, t,u) = -(l- 2e) 2u 2 - e(t 2 + s 2 + u 2 ) + 3sV + sV Box b (s, u) 
s 

\Qu 2 - 2t 2 e - e 2 {2t 2 + 5u 2 + 3tu) - eV 



ts(2 + e) 



+e 4 ts 



[Bub(w) - Bub(s)] . 



(4.26) 



Finally the square of the finite part of the one-loop amplitude is controlled by the 
finite functions and f± 



1 



f 3 (s,t,u) = - {(s 2 + w 2 )Box 6 (s,t) + (2m- s)[Bub(s) -Bub(t)]} + (9(e), (4.27) 



and 



h(s,t,u) = ^{2sBox e (t,u) + 3[Bub(w) -Bub(t)]} + 0(e). (4.28) 

Again, the infrared singular structure obtained by explicit expansion of T1Z t and 
XlZ nt as series in e, agrees with that obtained using Catani's formalism [13] 



on 



t,c 



r(l-e) Ve 2 2eJ [N 



"7 + ^ T7 ~ +TT -— " 



(N 2 + 1) / 



2\ e 



N 



and 



3^7 



nt,C 



T(l - e) Ve 2 2e 



TV 2 - 1 
N 



s 



1 

iV 



N 



+ 



u 



(4.29) 



(4.30) 

As before, we can rewrite Eq. (4.18) directly in terms of 11Z t ,c and lTZ n t,c rather 
than XTZ t and TlZ t ,c provided the finite difference is evaluated through to 0(e 2 ). 
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5. Conclusions 



In this paper we discussed the 0(af) virtual corrections to like and unlike mass- 
less quark-quark scattering formed by the interference of one-loop amplitudes with 
one-loop amplitudes. Our main results are Eqs. (4.1) and (4.18) where we provided 
analytic formulae for the MS-renormalised amplitudes in terms of two one-loop mas- 
ter integrals, the box graph in 6 — 2e dimensions and the bubble diagram in 4 — 2e 
dimensions. Expressions for the crossed and time reversed processes can be simply 
produced by inserting the expansions of the master integrals in the appropriate kine- 
matic region. Together with the interference of two-loop diagrams with tree graphs 
computed in [11, 12], the one-loop square contributions given in Eqs. (4.1) and (4.18) 
complete the set of 2 — > 2 virtual corrections for massless quark-quark scattering at 
NNLO. 
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